Analysis of Cirrus Cloud Microphysical Data 


Final Report 
NASA Grant NAG1-1706 


Michael R. Poellot 
Cedric A Grainger 


Department of Atmospheric Sciences 
University of North Dakota 
P.O. Box 9006 

Grand Forks, ND 58202-9006 


October 15, 1999 



TABLE OF CONTENTS 


Page 

1.0 Introduction 1 

2.0 Microphysical Data 2 

3.0 Approach 

3.1 Microphysical Characterization of Cirrus 3 

3.2 Analysis of the Spatial Variation of Cirrus Clouds 4 

4.0 Results 

4. 1 Cloud Microphysics 5 

4.2 Spatial Analysis 5 

References 1 3 


Appendix 


15 



1. Introduction 


The FIRE program has the goal of improving our capabilities to understand, 
model and detect the properties of climatically-important clouds. This is being 
undertaken through a three-pronged effort of modeling, long-term observations and short- 
term intensive field studies. Through examination of satellite and other data it is apparent 
that stratus and cirrus cloud types have the greatest impact on climate due to their 
radiative effects and ubiquitous nature. As a result, the FIRE program has developed two 
paths of investigation, each having its own subset of research objectives and 
measurement programs. The work conducted under this grant was directed toward 
furthering our understanding of cirrus cloud systems. 

While it is known that cirrus are climatically important, the magnitude and even 
sign of the impact is unclear. Cirrus clouds affect the transfer of radiation according to 
their physical depth and location in the atmosphere and their microphysical composition 
(Liou, 1986). However, significant uncertainties still exist in how cirrus clouds form and 
how they are maintained, what their physical properties are and how they can be 
parameterized in numerical models. Better remote sensing techniques for monitoring 
cirrus cloud systems and improved modeling of radiative transfer through ice particles are 
also needed. A critical element in resolving these issues is a better understanding of 
cirrus cloud microphysical properties and how they vary. 

Important microphysical parameters include ice water content and integrated ice 
water path, particle size distribution and crystal habit; the significance of these 
parameters has been shown in a number of applications. For example, Starr and Cox 
(1985) found notable interactions between cloud processes and microphysical properties 
in their modeling studies. They concluded that the bulk physical properties of cirrus are 
affected by crystal habit and size distribution, and that cirrus cannot be represented by 
analogy to warmer stratiform clouds. Stackhouse and Stephens (1994) illustrated how 
biases in particle size parameterizations can lead to biases in model results and 
interpretations. Model results of Mitchell and Amott (1994) showed a strong dependence 
of cirrus cloud radiative properties on particle habit. Significant differences resulted 
from the use of non-spherical particles versus equivalent spheres. Another critical area is 
the retrieval of cloud properties from remote measurements which requires some 
assumptions concerning the composition of the clouds. For example, Intrieri, et al, 

(1993), have developed a method for determining cirrus cloud particle sizes using lidar 
and radar measurements based on a modified gamma distribution of particle sizes. Other 
retrievals such as those using satellite data are also dependent on similar assumptions 
(e.g., Ackerman, et al, 1993). 

Our understanding of cirrus cloud microphysical characteristics has lagged to a 
large extent because of the relative inaccessibility of these upper-tropospheric cloud 
systems and a lack of adequate instrumentation for their measurement. Prior to 1991, 
only a limited number of in situ measurement programs had been conducted, and these 
included a variety of cirrus types and geographical locations. Perhaps not surprisingly, 
analyses of these data sets have found a high degree of variability in particle sizes and 
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concentrations (for reviews, see Liou, 1986, and Dowling and Radke, 1990). While there 
have been some parameterizations derived for modeling use (Heymsfield and Platt, 1984, 
Stephens, et al, 1990), there remains a significant amount of uncertainty as to the source 
and nature of this variability. This is particularly true at the small end of the size 
spectrum because most of the particle probes used in the past could not detect particles 
smaller than approximately 50 pm. There is evidence which suggests that there are 
enough ice particles in cirrus at sizes smaller than this to have a significant radiative 
impact (e.g., Platt, et al, 1989). The shape of the particle spectrum in this range is also 
needed for accurate retrievals of cirrus properties from remote sensors. Fortunately, the 
microphysical data collected in situ during the FIRE Cirrus IFO II included 
measurements of small ice particles. Analysis of these data was conducted to help reduce 
the uncertainties in cirrus microphysical properties. 

The focus of the research to be conducted under this grant was the data collected in 
situ by the University of North Dakota Citation aircraft. The goals of this research were 
to: 


• add to the body of knowledge of cirrus cloud microphysics, particularly at the small 
end of the size spectrum; and 
•analyze the spatial variation of cirrus clouds. 


2. Microphysical Data 

The FIRE Cirrus IFO II was conducted over southeastern Kansas in late fall of 
1991 , to sample clouds of a more sub-tropical origin and to expand upon the knowledge 
gained from the first IFO. Sampling was conducted simultaneously by active and passive 
remote sensors and by in situ aircraft and sonde-bom instruments. The mission of the 
UND Citation aircraft was to collect cloud microphysical and supporting data. A 
compliment of cloud particle sizing probes included a continuous formvar replicator and 
four standard Particle Measuring Systems (PMS) instruments: 


FSSP 

Cloud water droplets 

3 - 47 pm 

ID-C 

Cloud particles 

20 - 600 pm 

2D-C 

Cloud particles 

33 - 1056 pm 

1D-P 

Large particles 

0.3 - 4.5 mm. 


The replicator, supplied by the University of Nevada Desert Research Institute (DRI), is 
capable of capturing particles as small as approximately 7 pm. The FSSP is a light 
scattering instrument and its response to ice particles is uncertain (but it still may be 
possible to use its data to derive small particle information under some conditions). The 
other three are optical array probes. Liquid and supercooled water contents were 
measured with a Johnson-Williams probe and Rosemount ice detector, respectively. 

Other measured parameters included the state parameters (temperature, pressure and dew 
point), three-dimensional wind flow and turbulence, condensation nuclei (CN), ozone and 
NO x concentrations, and aircraft dynamics. 
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A variety of mid and upper-tropospheric cirriform clouds were sampled during 
the IFO. A thorough overview of the conditions encountered is given in a technical 
report by Poellot and Henderson (1994); the highlights may be summarized as follows. 
Ten missions were flown for a total of 35 flight hours at cloud altitudes ranging from 
approximately 6.6 to 12.2 km with in-cloud temperatures of - 23 °C to -64 °C. The 
typical flight profile consisted of step climbs and descents on an along-wind track; 
several spiral descents were also performed. Cirrus conditions sampled included layers 
which were optically thin, optically thick, "patchy", multi-layered, deepening or 
dissipating. Supercooled water was encountered on two flights. Many of the clouds 
appeared to contain an abundance of small (<100 pm diameter) ice particles. 

A recent analysis of a portions of these data compared the performance of the 2D- 
C optical array probe, which is the type most often used in previous measurement 
programs, and the formvar replicator (Amott, et al, 1994). The comparison of 2D-C and 
replicator data revealed that for this sample cloud there were some tens to hundreds of 
small particles which went undetected by the 2D-C. Calculations showed that these 
particles could significantly contribute to and even dominate the solar extinction and 
infrared emission of cirrus. Therefore, the replicator data were used in the analyses to 
extend our knowledge of cirrus microphysics to the smaller sizes. 


3. Approach 

3.1 Microphysical Characterization of Cirrus 

Ice water content, particle size distribution and habit are microphysical 
parameters important to the proper modeling and monitoring of cirrus clouds and their 
radiative effects. The method of analysis followed that of Heymsfield and Platt (1984; 
hereafter referred to as HP), where pass average values of particle concentrations by size 
were calculated. Data from the DRI replicator were used to count and size particles with 
diameters up to 150 pm and the 2D-C for larger ones. Crystal habits were determined 
from the replicator data for particles smaller than 300 pm and from the 2D-C data for 
larger ones through the use of an automatic classification scheme (Heymsfield and 
Pamsh 1979). (Particles larger than about 300 pm tend to fragment when impacting the 
replicator film surface.) Ice water contents (IWC) were calculated using the pass- 
averaged concentrations and crystal habits using relationships cited in Heymsfield (1977) 
and Mitchell (1994). 

Following HP, the size distributions were parameterized for each 5 °C 
temperature range analyzed, and the pass average IWC used to normalize the spectra. 
Values obtained for the clouds sampled in this study were compared with those derived 
by HP. Some differences were expected due to the cloud types sampled and the 
instrumentation used to collect the data. The data set used by HP was taken from 
primarily deep, uniform cirrus generally associated with synoptic depressions, while the 
FIRE samples were mostly thin, banded features associated with jet stream circulations or 
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baroclinic ascent. Also in contrast with the present study, the small particle size spectra 
(down to 20 pm) were derived by HP from ID probe data. 


In addition to the stratification by temperature, an effort was made to determine 
relationships between the microphysical and macrophysical cloud properties such as 
cloud depth, distance from cloud boundary, and large-scale forcing. A review of 
complimentary aircraft, lidar and radar data was performed to determine the study 
clouds’ top and base heights and the relative position of each Citation sample leg within 
cloud. Supporting synoptic descriptions and satellite data were also reviewed to 
determine the large scale forcing of the sampled cirrus cloud systems. 


Analysis of the replicator data for this study was accomplished under subcontract 
with DRI. The tedious job of counting and sizing particle replica has been somewhat 
streamlined by the acquisition of equipment and development of software to make the 
process semi-automatic (Amott and Hallett, 1994). Even so, it is still a time-consuming 
procedure and on the order of a half million frames of replicator film were exposed 
during the FERE EFO alone. Therefore, it was not within the scope of this proposal to 
analyze all of the replicator data; rather, portions of the film from a representative portion 
of each data leg were processed. 


3.2 Analysis of the Spatial Variation of Cirrus Clouds. 

The variation of the microphysical characteristics of the cirrus clouds over 
distances of a few meters to tens of kilometers is important when trying to interpret data 
with relatively coarse resolution, such as radiometric measurements of cloud 
characteristics. The manner in which the cirrus microphysical characteristics vary would 
be expected to be closely related to the dynamics of the generating mechanism. Sassen et 
al. (1989) have shown by a spectral analysis of aircraft-derived vertical velocities that 
there is an apparent convective structure present in cirrus with a horizontal scale of 
approximately 1-2 km. Therefore, an analysis of the vertical velocities derived from the 
UND Citation data for the legs flown in the FIRE cirrus was accomplished. The purpose 
of this analysis was to see if the type of convective structures reported by Sassen et al. are 
apparent in this larger data base. In addition, similar analyses were performed on several 
of the cirrus microphysical parameters. The parameters included ice crystal 
concentration, crystal size, and ice water content. The output of this work was a 
description of the variation of these parameters as a function of scale, to see if the 
apparent fluctuations in vertical velocity at the 1 km scale are indicative of fluctuations in 
cloud microphysical characteristics of the same scale. 

In a similar vein, if there are sustained convective structures responsible for the 
peaks in the spectral density described earlier (rather than random fluctuations on that 
scale) there would be a reasonable likelihood of finding a significant positive eddy flux of 
ice crystals (or of ice water content) that is independent of the larger scale (i.e., >100 km) 
vertical motions. This would be expected since a sustained downdraft would contain 
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crystals that are sublimating while the updraft would contain crystals that are growing. 
Non-sustained vertical motions, lasting for relatively short time periods would result in 
very little vertical displacement and any effects would not be likely to be seen in the 
microphysics data. Therefore, the eddy flux computations should be near zero. The 
presence or absence of sustained convective elements within the cirrus would be of major 
importance in understanding the microphysical characteristics of these clouds 


4. Results 

4.1 Cloud Microphysics 

The detailed results of the cirrus cloud microphysics are given in Poellot, et al., 
1998. The following is a summary of these results. 

This analysis of the cloud microphysical data collected during FIRE II by the 
UND Citation found that there are a number of factors which influence cirrus cloud 
microphysics and which must be taken into account for parameterization of these cloud 
systems. These include but are not necessarily limited to temperature, synoptic forcing, 
position in cloud and cloud ice water content. The key results are summarized as 
follows: 1) A dependence of particle size spectra shape on temperature was found which 
differed from that previously reported for deeper, frontal overrunning cirrus. 2) There 
was no marked dependence of crystal habit on temperature other than columns were most 
prevalent in clouds colder than - 50°C. 3) When synoptic forcing was considered, 
clouds generated by closed low systems and the sub-tropical jet (southwesterly flow 
cases) contained many more particles at all sizes than those produced by short wave and 
baroclinic lift (northwesterly flow cases). 4) A position-in-cloud stratification supported 
earlier studies with evidence that ice particles are nucleated in the upper portion of the 
cirrus cloud, grow by diffusion and aggregation as they fall, and finally evaporate as they 
reach cloud base. 5) A look at the variation of spectral form with ice water content 
(IWC, a common model-predicted variable) found that as IWC increased, an increase in 
particle number concentration was accompanied by a broadening of the distribution to 
larger sizes. 6) Analysis of the replicator data suggests that the small particle end of the 
size spectrum may be best represented by a modified Gamma distribution. 

4.2 Spectral Analysis 

The spectral analysis performed on each of the components of the wind velocity 
showed that, with few exceptions, the shape of the spectra was similar and the turbulence 
intensities were comparable. The slopes of the power spectra for each of the wind 

velocity components were close to -j for the high wavenumber end of the spectrum. The 

vertical velocity spectra were run on most of the days. Some examples of the vertical 
component spectra are shown in figures 1-3. These spectra are unfiltered, so there are no 
known biases for any particular portion of the spectrum. On many of the time blocs 
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investigated, the spectra showed a significant change in the slope for wavenumbers 
corresponding to wavelengths of between 1 and 2 km (wavenumbers of between .003 and 
.006 m' 1 ). Generally, the slope became closer to zero at the very low wavenumbers. This 
would indicate that there was an input of energy at the 1-2 km scale. 

The spectral analyses on the microphysics data showed little consistent difference 
from one parameter to another in terms of the shapes of the curves. That is, the power 
spectra for the 2D-C concentrations were not markedly different from the spectra of the 
1D-C concentrations. Some typical spectral density curves for the 1D-C probe are shown 
in figures 4-6. These are for data taken in the same time periods as was done for the 
vertical velocity spectra in figures 1-3. In general, the spectra of the hydrometeor 
concentrations tended to be fairly flat (i.e., near zero slope) at the high wavenumber end 
of the spectrum. This “whiteness” of the spectrum implies that there are no preferences 
for variations in the microphysical data for one size scale to another. There was 
substantial variability at the low wavenumber end of the spectrum and some variability in 
the point where the slope of the spectra began to deviate from zero. There was also 
significant day to day differences in the magnitude of the variations as indicated by the 
spectral magnitudes. 
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Figure 1 . Spectral density for the vertical velocity as a function of wavenumber for the 
period 16:10:00 to 16:23:50 on 14 November 1991. The dashed lines represent a -5/3 
si 



0 . 1 0E+00 


0. 10E-01 


0. 10E-02 


0.10E-03 


0.10E-04 


0.10E-05 


0.10E-06 

0 


511 


Power Spectrum Density 


1 ^ 1 1 riTT| I ] 1 TTTT] i ^1 1 I 1 ITT] I T" 



0001 


0.0010 0.0100 0.1000 

Wave Number (Per Meter) 

11 / 22/91 


Figure 2. Same as Figure 1, but for the time period 18:30:40 to 18:34:20 on 22 November 
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Figure 4. Spectral density for the 1D-C concentrations as a function of wavenumber for 
the period 16:10:00 to 16:23:50 on 14 November 1991. The dashed lines represent a -5/3 
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Figure 5. Same as Figure 4, but for the time period 18:30:40 to 18:34:20 on 22 November 
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Figure 6. Same as Figure 4, but for the time period 16:27:00 to 16:30:00 on 30 November 
1991. 
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Appendix: Particle Size Spectra 

This appendix contains a series of graphs showing detailed particle size spectra 
and distribution parameters derived from the Citation microphysical data. The data 
consisted of aircraft pass-averaged spectra, which were divided into temperature 
categories and location in cloud categories. The temperature categories were intervals of 
5 C, the location categories were either bottom, bottom-middle, middle, top-middle, or 
top. These data are summarized in the attached figures. 

In the case of the 1DC and 2 DC data set, exponential distributions were fit to the 
data in a manner similar to that used by Heymsfield and Platt (1984). The parameters 
they discussed are compared with the data collected here and summarized in the attached 
figures and table. 

The empirical corrections made to the 1DC data were then removed and a gamma 
distribution was fit to this uncorrected 1DC and 2DC data. These can be found in the 
attached figures for both temperature and location in cloud categories. The figures titled 
“1DC Gamma Distribution Parameters vs. Temperature” and “1DC Gamma Distribution 
Parameters vs. Location in Cloud” give a graphical representation of the parameters 
obtained from the best-fit gamma distribution. The parameter labeled Gamma is 
representative of the ‘average’ concentration in the distribution and the parameter labeled 
Peak is representative of the size at which the maximum concentration in the distribution 
occurred. Similar figures are given for gamma distributions fit to the replicator and 2DC 
data. 


Replicator ice particle habits: the ice particle habit data from the replicator was 
also obtained and partially analyzed. Two figures are included from this limited analysis. 
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Particle Size Distributions 


by Temperature and Location in Cloud 



Temperature Averaged Spectra, Replicator & 2DC Probe 

(-35°C to -40°C, N = 7) 



o 

r— ( 

<N 

co 



VO 

c- 

o 

O 

o 

O 

O 

o 

o 

o 

+ 

W 

w 

tq 

UJ 

W 

W 

W 

tq 


1 

(——I 


i— • 

— 

— h 



((uojoiui J3JH)/#) ap/NP 


Diameter (microns) 




Temperature Averaged Spectra, Replicator & 2DC Probe 

(-40°C to -45°C, N = 4) 



o 

f— i 

rN 


Tj- 

in 

VO 

r- 

o 

o 

o 

o 

o 

o 

o 

o 

“h 

W 


W 

W 

W 

tq 

W 

i 

» — 1 

— • 





- — ' 


((UOJDIUI QP/KP 


Diameter (microns) 




Temperature Averaged Spectra, Replicator & 2DC Probe 

(-45°C to -50°C, N = 6) 



o 

r— ' 


m 

rf 



r- 

o 

o 

o 

o 

o 

o 

o 

o 

+ 

W 


Li 

W 

ti 

ti 

ti 

Li 



— 

i— • 




— 


((UOJOIUI J9JII)/#) OP/NP 


Diameter (microns) 




Temperature Averaged Spectra, Replicator & 2DC Probe 

(-50°C to -55°C,N= 15) 



Diameter (microns) 




Temperature Averaged Spectra, Replicator & 2DC Probe 

(-55°C to -60°C, N = 6) 



1— > 

CN 

on 



'O 

r- 

o 

1 

O 

O 

o 

o 

o 

o 

tq 

oq 

w 

ti 

W 

PJ 

td 

t— 

— 







((uOJOtUI J3U[)/#) QP/NP 


Diameter (microns) 




Temperature Averaged Spectra, Replicator & 2DC Probe 

(-60°C to -65°C, N = 2) 



o 

i— < 

cm 

cc 

^r 


vO 

r- 

o 

o 

o 

o 

o 

o 

o 

o 

“h 

W 

UJ 

tu 

U 

PJ 

ti 

ci 

ci 



1 


* 

•—* 

— 

■— < 


((UOJOIUI J91l[)/#) OP/KP 


Diameter (microns) 













Exponential Distribution Fit to Size Spectra 





00 

05 

T— 

oo 

CO 

CO 


M- 

M- 



CD 

00 

m 

CO 

oo 

o 

m 


I s - 



CM 

o 

in 

CO 


■*— 

co 

q 

tT 



CO 

i 

CO 


in 

d 

i 

CO 

00 

1 

r-^ 

i 

od 

1 

CM 











GQ 


CO 


M" 

CO 

05 

CO 

CO 

co 

in 

o 


m 

t— 

OO 

T— 

CO 

o 

CO 

05 

o 

S 

o> 

I s - 

in 

co 

05 


o 

CM 

i 

UJ 


*r 

M- 

in 

m 

o 

O 

05 

oo 

O 


E 

oo 

I s - 

00 

CO 

O 

I s - 

CM 

I s - 

r— 

o 

n 

o> 

05 

o 

o 

CO 

o 

o 

O 

CO 

in 

tn 



o 

o 

o 

O 


Z 

3 

in 



o 

d 


d 

o 

05 



m 

00 

CM 

CO 

r- 

CM 

m 

m 

oo 


*7 

CM 

o 

CM 

I s - 

o 

oo 

q 

o 

o 


E 

’T 

v- 

CM 

q 

CO 

o 

o 

M* 

o 

UJ 

UJ 

1 

UJ 

Q 

3. 

o 

O 


o 

o 

o 

r- 

CO 

oo 

g 

r> 

g 



o 

d 

o 

q 

r— ; 

q 

q 

z 






d 


CO 

r-^ 



CO 

05 

T— 

CO 

co 

r- 

CM 

CO 

CM 




CO 

T~ 

T— 

m 

00 

t — 

oo 

Tf 



05 

q 

q 

q 

o 

00 

q 

oo 




T— 

r— 

T— 


05 


d 

d 

d 



1 

i 

1 

1 

d 

d 

i 

i 

i 

-r- 











GQ 


m 

CO 

m 

r- 

in 

m 

05 

I s - 

T— 

O 


M* 

r 

m 

I s - 

•r- 

r— 

t— 

CM 

o 

T ^~ s 

00 

05 

oo 

05 

CO 

OO 

I s - 

o 

CO 

q 

I 

E 

M- 

Tf 

o 

05 

d 

00 

05 

od 



d 

in 

d 

1^ 

m 


05 

oo 

co 


— 

CO 

I s - 

CM 

M- 

M' 

I s - 

05 


I s - 

o 

T" 

05 

m 

CM 

CO 

m 

r— 

in 

to 

oo 

Z 

3 

CM 

CM 

M- 

T— 

yr~ 


T ” 





T— 

m 

CO 

M; 

T— 

m 

T- 

q 

CO 



CO 

d 

CO 

M 1 

in 

CD 

o 

co* 

d 


E 

CO 

CO 

00 

CM 

I s - 

o 


oo 

M" 





CM 

CM 

t— 















o 

o 










z 

E, 

c- 

T— 

in 

r- 

m 

CO 

CO 

M* 

r- 

o 


CM 

05 

M* 

M- 

o 

05 

T — 

CO 

CM 

1 


CO 

tT 

o 

CO 

00 

o 

CO 

T— 

OO 

op 

CM 

CO 

CM 

co 

I s - 

CO 

co 

m 

o 


CM 

CM 

o 

T— 

T— 

o 

o 

o 

o 

O) 

E 

O 

O 

d 

o 

o 

q 

o 

o 

o 

> 

< 

3 

o 

o 


d 

o 

d 

o 

o 

o 



r-- 

00 

oo 

I s - 

OO 

m 

co 

M* 

CO 

o 


05 

M- 

oo 

o 

OO 

CO 

o 

CM 

00 


r- 

00 

CM 

m 

m 

o 

CM 

m 

T— 


cp"^ 

o 

t— 

o 

o 

o 

o 

o 

o 

o 



T~ 

o 

o 

o 

q 

o 

o 

o 

c 

E 

O 

o 

o 

o 

o 

d 

o 

o 

q 

5 

3 

o 

d 

o 

o 

o 


o 

o 

d 

O 


O) 

05 

CM 

m 

in 

m 

m 

m 

00 


M* 

M- 

^ — 

CO 

oo 

T“ 

M* 

co 

co 

§ 


r- 

m 

M- 

CD 

CO 

o 

CM 

r-- 

M* 

CO 

05 

CO 

CO 

CO 

o 

in 

CO 

CM 

— 


CO 

CO 

oo 

00 

in 

T— 

CM 

CM 

o 

X 

CO 

5 

E 

O 

o 

o 

o 

o 

o 

O 

o 

o 

3 

o 

d 

d 

d 

o 

o 

o 

o 

d 

{/) 


M- 

m 

CO 


CM 

CM 


oo 

CM 

<D 





CM 

T— 


CM 



a. 











E 











CO 











CO 











^ , 











o 











t 











0) 











.a 











E 











3 











z 











a> 











o> 











c 


to 

o 

m 

O 

in 

o 

in 

o 

in 

CO 

cr 

Q_ 


CM 

CO 

< 

CO 

i 

TT 

i 

M- 

i 

m 

i 

m 

1 

CO 

1 

co 

1 


O 

o 

o 

o 

O 

o 

o 

o 

o 












E 


o 

m 

o 

m 

o 

in 

o 

in 

o 

(D 


CM 

CM 

CO 

co 

M- 


m 

m 

CO 

1 

H 


i 

• 

1 

i 

i 

i 

i 

i 



Temperature Averaged Spectra, 1DC & 2DC Probes 

(-20°C to -25°C, N = 4) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-30°C to -35°C,N= 16) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-35°C to -40°C, N = 21) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-40°C to -45°C,N= 12) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-45°C to -50°C,N= 12) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-50°C to -55°C, N = 21) 
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Temperature Averaged Spectra, 1DC & 2DC Probes 

(-55°C to -60°C, N = 8) 
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Temperature Averaged Spectra, 1DC & 2D C Probes 

(-60°C to -65°C, N = 2) 
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(T1 = -20°C, DT = 5°C) 
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Gamma Distribution Fit to Size Spectra - 1DC and 2DC 
by Temperature 






eranrer) 





Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

. (-20°C to -25°C, N =4) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-25°C to -30°C, N =5) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-30°C to -35°C, N =16) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-35°C to -40°C, N = 21) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-40°C to -45°C,N= 12) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-45°C to -50°C,N= 12) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best- 

Distribution 

(-50°C to -55°C, N = 21) 
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Temperature Averaged Spectra (1DC & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-60°C to -65°C, N = 2) 



Diameter (microns) 






Gamma Distribution Fit to Size Spectra - 1DC and 2DC 
by Location in Cloud 
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Gamma Distribution Fit to Size Spectra - Replicator and 2DC 


by Temperature 



Replicator Gamma Distribution Parameters vs Temperature 
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Temperature Averaged Spectra (Replicator & 2DC Probe) and Best-Fit Gamma 

Distribution 
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Temperature Averaged Spectra (Replicator & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-40°C to -45°C, N = 5) 
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Temperature Averaged Spectra (Replicator & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-45°C to -50°C, N = 6) 
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Temperature Averaged Spectra (Replicator & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-55°C to -60°C, N = 6) 
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Temperature Averaged Spectra (Replicator & 2DC Probe) and Best-Fit Gamma 

Distribution 

(-60°C to -65°C, N = 2) 
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Gamma Distribution Fit to Size Spectra - Replicator and 2DC 
by Location in Cloud 
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Location in Cloud Averaged Spectra (Replicator & 2DC Probe) and Best-Fit 

Gamma Distribution 

(bottom, N = 6) 
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Gamma Distribution 
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Location in Cloud Averaged Spectra (Replicator & 2DC Probe) and Best-Fit 

Gamma Distribution 

(top-middle, N = 7) 
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Location in Cloud Averaged Spectra (Replicator & 2DC Probe) and Best-Fit 

Gamma Distribution 

(top, N = 9) 
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Crystal Habit Information 
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